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Abstract 

A further development of the evolutionary picture of A+A collisions, which we call the inte¬ 
grated HydroKinetic Model (IHKM), is proposed. The model comprises a generator of the initial 
state GLISSANDO, pre-thermal dynamics of A+A collisions leading to thermalization, subse¬ 
quent relativistic viscous hydrodynamic expansion of quark-gluon and hadron medium (vHLLE), 
its particlization, and hnally hadronic cascade ultrarelativistic QMD. We calculate mid-rapidity 
charged-particle multiplicities, pion, kaon, and antiproton spectra, charged-particle elliptic flows, 
and pion interferometry radii for Pb-|-Pb collisions at the energies available at the CERN Large 
Hadron Collider, ^/s = 2.76 TeV, at different centralities. We find that the best description of the 
experimental data is reached when the initial states are attributed to the very small initial time 
0.1 fm/c, the pre-thermal stage (thermalization process) lasts at least until 1 fm/c, and the shear 
viscosity at the hydrodynamic stage of the matter evolution has its minimal value, g/s = ^. At 
the same time it is observed that the various momentum anisotropies of the initial states, different 
initial and relaxation times, as well as even a treatment of the pre-thermal stage within just viscous 
or ideal hydrodynamic approach, leads sometimes to worse but nevertheless similar results, if the 
normalization of maximal initial energy density in most central events is adjusted to reproduce 
the final hadron multiplicity in each scenario. This can explain a good enough data description in 
numerous variants of hybrid models without a prethermal stage when the initial energy densities 
are defined up to a common factor. 

PACS numbers: 25.75.-q, 24.10.Nz 
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I. INTRODUCTION 


Hydrodynamics is considered now as the basic part of a spatiotemporal picture of the 
matter evolution in the processes of ultrarelativistic heavy ion collisions (see recent reviews 
l|, l2|). To complete the description of A+A collision processes, hydrodynamics must be 
supplied with a generator of an initial non-equilibrated state, pre-thermal dynamics which 
forms the initial near locally equilibrated conditions for hydro-evolution, and prescription 
for particle production during the breakup of the continuous medium at the final stage of 
the matter expansion. 

As for the initial state, since it fluctuates on an event-by-event basis, Monte Carlo event 
generators are widely used to simulate it in relativistic A + A collisions. The most com¬ 


monly used models of initial state are the 
kLN (Monte Carlo Kharzeev-Levin-Nardi) 


C-Glauber (Monte Carlo Glauber) [^, MC- 
EPOS (parton-based Gribov-Regge model) 
EKRT (perturbative QCD -|- saturation model) jb, 7|, and IP-Glasma (impact parameter 
dependent glasma) 8(]|j The last model also includes non-trivial non-equilibrium dynamics 
of the gluon fields which, however, does not lead to a proper equilibration. To apply these 
models to data description some thermalization process has to be assumed. Evidently, in 
order to reduce uncertainties of results obtained by means of hydrodynamical models, one 
needs to convert a far-from-equilibrium initial state of matter in a nucleus-nucleus collision 
to a close to locally equilibrated one by means of a reasonable pre-equilibrium dynamics. A 
relaxation time method 10|, initially developed for the post-hydrodynamic stage, is modified 
for applications to pre-thermal dynamics. 

As for the post-hydrodynamic (afterburner) stage, good progress was made in under¬ 
standing and modeling the breakup of hadron matter into hadron gas at the late stage of 
expansion. Typically, particlization, i.e., the transition from the hadron (or quark-gluon) 
fluid to hadronic gas expansions, is described by means of the so-called Cooper-Frye pre¬ 
scription within hybrid models. In such an approach, the sudden conversion of the fluid 
to particles happens at a hypersurface of hadronization or chemical freeze-out. It has long 
been known that such a matching prescription has problems with the energy-momentum 
conservation laws when fluid is converted to particles at the hypersurface which contains 


non-space-like parts (see up-to-date review [l| and references therein). These problems can 


^ Note that the last two models are able to reproduce correctly the centrality systematics of the Sv 2 7|, l9|. 
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be avoided by using the HydroKinetic Model (HKM) that was proposed in Ref, 


further developed in Refs. 




n. 


Cl 


and 


(see also Ref. [l^). In basic HKM, particlization is 
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In 


considered not as a sudden, as in hybrid models, but as a realistic continuous process 
such an approach the description of continuous particle emission from the fluid is based not 
on the distribution function but on the so-called escape function InlQ. 

In Refs. Ifl 16| the method of eseape probabilities, similar to that used in HKM for- 
malism for post-thermal dynamics, is applied to the pre-thermal one. It makes it possible 
to describe the pre-thermal evolution of the energy-momentum tensor using three free pa¬ 
rameters: the time of the initial (non-equilibrium) state formation, mean relaxation time, 
and time of thermalization, say, 1 fm/c. This phenomenological model allows one to provide 
initialization of the hydrodynamical evolution, and is a more general and simple alternative 
to the so-called “anisotropic hydrodynamics” 17 


2l|, as discussed in detail in 1^ . 


The target energy-momentum tensor of the pre-thermal dynamics of A+A collisions, 
which the initially created system reaches during evolution, is the tensor of viscous relativistic 
hydrodynamics. The further evolution is governed by the equations of the Israel-Stewart 
framework of relativistic viscous hydrodynamics until the lowest possible temperature when 
the system is still close to local thermal and chemical equilibrium. It dehnes the hypersurface 
which is often called the chemical freeze-out isotherm. The particlization procedure starts 
typically at this very hypersurface. The hypersurface contains non-space-like parts that 
formally enforce one to use the HKM continuous particlization procedure, to eliminate the 
problem with violation of the energy-momentum conservation. However, with almost smooth 
(averaged) initial conditions, the effects originated from these parts are quite small and do 
not affect signihcantly the observables 13(]lj Moreover, in viscous hydrodynamics the effects 
of violation of local equilibrium, which is produced in HKM for an ideal fluid, are partially 
taken into account. Therefore, summarizing all that and accounting for very time consuming 
calculations in HKM, we will use in this paper the sudden transition from hadron matter to 
hadron gas, similar to the hybrid models. In future studies we plan to include the option 
of continuous HKM particlization for situations when undesirable contributions from non¬ 
space-like parts are fairly large: it can be for dvr- geometry observables and event-by-event 


analysis. Considering also that an important component of the 


nil model is the pre-thermal 


^ For recent discussions of the particlization procedure see, e.g., Refs. [l4 IlSjk 
^ The reason is that at small proper times r the three-dimensional areas of non-space-like parts are small, 


and for large times the hydro-velocities at the corresponding parts are large. 
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evolution where the description is built on the same idea as the basic HKM formalism, we 
will call the complete model, which includes generation of the initial state, its thermalization, 
sudden (with the possible option of “continuous”) particlization, subsequent ultrarelativistic 
QMD (UrQMD) hadronic cascade 2^, and one- and multi- particle spectra formation, the 
integrated HydroKinetic Model (iHKM). 

The paper is organized as follows. Section 2 includes subsections that describe the ingre¬ 
dients of the complete iHKM. Section 3 is devoted to the results and discussions, while the 
Section 4 summarizes the physical meaning of obtained results. 


II. THE MODEL DESCRIPTION 

As discussed in the Introduction the integrated HydroKinetic Model (iHKM) consists of 
hve ingredients which we describe below. 


Initial state 


Though prethermal dynamics is used to initialize hydrodynamics, it itself needs the initial 
state of the matter at the starting time tq, say tq = 0.1 or 0.5 fm/c, when one can speak about 
energy density distribution in non equilibrated matter. Following Ref. 23|, an appropriate 
scale controlling the formation of gluons with a physically well-defined energy is roughly 
3/(7^/i, where = Airag and /i^ are dimensionless parameters, which is the variance of the 
Gaussian weight over the color charges p of partons. The estimate of this time gives Tq = 0.1 
fm/c for the top energy available at the BNL Relativistic heavy Ion Collider (RHIC) and 
To = 0.07 fm/c for the CERN Large Hadron Collider (LHC) full energy 24 1. We will utilize 
To = 0.1 fm/c as the default value for the LHC energy ^ysNN = 2.76 TeV; however, we will 
use also Tq = 0.5 fm/c to analyze the corresponding variation of the results. 

We attribute the initial state to the time tq. The generation of the initial state is based 
on the GLISSANDO 2 25|, l26| package. This package works in the frame of the semiclassical 


Glauber model. Within this approach at the very initial stage of the collision, individual 
interactions between nucleons deposit transverse energy. Each deposition of the transverse 
energy at a certain space-time point or region is called a source and each source has its 
weight which is called relative deposited strength, or RDS. The normalization of RDS can 
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be treated as an insertion of the additional parameter. We choose its valne in such a way 
that all charged-particle multiplicity distribution, obtained in iHKM for 5% most central 
collisions, hts its experimental value. The RDS may be different for the wounded nucleons 
and the binary collisions, and it can fluctuate from source to source. We use a mixed model, 
amending the wounded nucleon model with some binary collisions. In this model, a wounded 
nucleon obtains the RDS of (1 — a)/2 and a binary collision has the RDS of a. The total 
RDS, averaged over events, is then (1 — a)iViy/2-|-aiVftm- Of course, the result of simulations 
fluctuates from event to event. In this paper we do not study fluctuations of observables, 
and will simulate inclusive observables basing on the approximation of the averaged initial 
state for each centrality class. This approximation implies that for each selected centrality 
only the one hydro run is performed. The simulation of a single event is made in three 
stages: 

1. Generation of the positions of nucleons in the two colliding nuclei according to the 
fluctuating nuclear density distribution. The form of the mean distribution depends 
on the mass of nuclei A. For sufficiently large nuclei this distribution has a Woods- 
Saxon form with an account of nuclear deformations, the latter are small in our case 
of Pb-I-Pb collisions. 

2. Generation of the transverse positions of the sources and their RDSs. 

3. Galculation of the physical quantities and writing the results in the output file. 

In order to obtain the initial conditions from the RDS, one can put the RDS to be propor¬ 
tional to energy or entropy. We choose it to be proportional to energy because entropy is not 
created yet in a non-equilibrium initial state. Then, for the initial state averaged over many 
events with a chosen centrality bin, one can write {b is the impact parameter associated with 
centrality) 

{1 - a)Nw{b,Tj,)/2 + aNhin{b,rT) . . 

et , To, rr) - = Q, = 0)/2 + aNUb = 0, = 0) ’ ^ ^ 

where eo is a fitting parameter dehning the maximal energy density of the initial state that 
is reached at minimal centrality, say 0-5% (5 ~ 0), at = 0. The value eo and a are the 
same for all centrality classes. The values of Cq and a are hxed by the reproduction of the 
charged-particle multiplicity at different centralities. 
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B. Pre-thermal stage of the matter evolution 


The generation of the initial energy density configuration with the help of GLISSANDO 
is supposed to be followed by a thermalization process. As we have already discussed in the 
Introduction we describe the matter evolution at the pre-thermal stage with the relaxation 
model 


10 


16| . We assume a longitudinal boost invariance, which is a good approximation 


for the central rapidity region of the fireball at high collision energies. The relaxation model 
simulations are applied to the system for the early stage of the heavy ion collision, from the 
proper time t = Tq to t = Tth, where r = y/t^ — As mentioned above, we use tq = 0.1 
fm/c as the basic value and tq = 0.5 fm/c for a comparison. Thermalization time is supposed 
to be Tth = 1-0 fm/c, the “conventional” time of the formation of strongly interaction matter. 


t ^ 1/m^EI The total energy-momentum tensor of the matter is taken in the form 
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16| 


( 2 ) 


where ^free(^) hydrodynamic (local equilibrium) and free (fully or almost 

free evolving) components of the energy-momentum tensor, V{t) is the weight function 
satisfying the following conditions: 0 < V{t) < 1, V{tq) = 1, Virth) = 0, and d^ViT)^^^ = 0. 
Its form will be discussed below. As for the free (or almost free) evolving component, it is 
well defined in such models as IP-Glasma and EPOS, but not in MG Glauber. To analyze 
different kinds of anisotropy of the initial state we present the boost-invariant distribution 
function at the initial hypersurface cxo: r = Tq in the factorized form 


/(Go, r^o, P) = e(^; T-o, rT)/o(p) (3) 

where e(5; tq, r^) is defined by Eq. ([T]). At other times /(x, p) is defined by the free-streaming 
requirement /(t, r, p) = /(Go,!" “ ^{t — Go),p), so that at all times the free evolution of 
the energy-momentum tensor is defined by the formula 

[d^p^f{x,p). ( 4 ) 

J Po 

Then the evolution of T/(/g(x) and corresponding distribution function f{x,p) are defined 
by the initial conditions, which are generated by GLISSANDO, and the initial momentum 
anisotropy of the function /o(p). For V{to) = 1, we have ^toL(D,x) = Tf((/;(ro,x). 

^ Note, that calculations with Tth = 1-5 fm/c and with the other parameters to be the same as in our basic 
scenario (see later) lead to very similar results as calculations with Tth = 1-0 fm/c. 
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The hydrodynamic component of the energy-momentnm tensor has the form 

^hydro(^) = [ehydro(x) + Phydro(a;) + ^Kydroi^Xydroi^) “ [Phydro(2:) + (5) 

where ehydro is the energy density in the local rest frame, Phydro is the pressnre, is the 
shear stress tensor, fl is the bnlk pressnre, and u’^ydmX i® energy flow fonr-vector. We 
neglect the bulk pressure, fl = 0. In curvilinear (hyperbolic) coordinates the equation of 
motion for the shear stress tensor is taken in the form 

-jrlJ-l' _ A 

= -^ ( 6 ) 

^TT ^ 

where the semicolon means the covariant derivation, brackets are defined as = 

-I- iA)(A^ — |A^^Aq,/ 3 )A"^, A'^'^ = — u^u'', and tins is the Navier-Stokes 

shear stress tensor: 


^Ns ^ p{Af^^d.\u'' + A'^^d-xX ~ IpA^^’^d-xu^ 


For V{t) we use the ansatz ^ 


(7) 


P(r) 


rth - T 


■^rel 


'Tth “ ^0 


( 8 ) 


where, as mentioned above, tq = 0.1 fm/c, Tth = 1-0 fm/c, and Trei < Tth — tq. The value 
of Trei is the model parameter, which characterizes the rate of thermalisation, that is, the 
speed of conversion of the non-equilibrium state to about the equilibrated one. Writing 
down the conservation laws for the total energy-momentum tensor in the form = 0 

and accounting that for the free streaming = 0, we have 


a^ll - nT)]T^^aJx)} = -TZ,(x)d,„V{,T). (9) 

Let us introduce the re-scaled hydrodynamic tensor with initial conditions Xy^^X) ~ 

0 at r = To for all x. Then 


= -ti‘;jx)b.,^v(t). ( lo ) 

This is the hydro dynamic-type equation of motion with a source. As mentioned above, 
7"/ree dehned by the the initial state at proper time Tq, and V{t) is dehned explicitly, 
then the source in the right-hand side of Eq. ca can be calculated for all r and x. The 
introduction of the re-scaled hydrodynamic tensor leads to the rescaling of the shear stress 
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tensor = 7r^'^(l — V). Then, multiplying Eq. ([6]) by (1 — V), we have the equation of 
motion for the re-scaled shear stress tensor 


[1 — P(r)] ( u^d- 


TT 


fll/ 


— [1 — V{t)]71\ 


(lU 

NS 


o 


( 11 ) 


’^a-nr))/ 

We also have to specify the equation of state in order to close the set of evolutionary 
equations. In all calculations the Laine-Schroeder equation of state (EOS) |27| is applied. 
At r = Tth we switch to the viscous hydrodynamic approach. The equations of motion are 
nearly the same as in Eq. m. but with the right-hand side equal to zero. 


C. Matter evolution in thermal and chemical locally near-equilibrated zone 


At T = Tth = ^ fm/c the V{Tth) = 0 and the target function is reached: = Tj^ydroi^)- 

The further evolution is described by the relativistic viscous hydrodynamics according to 
the equations ([5]) - ([7]), flTOl) . and flTT]) with 'P(r) = 0. A numerical solution of the viscous 


hydrodynamic equations is constructed with the vHLLE code 


Such an evolution describes 


the expansion of superdense quark-gluon and hadron matter close to local chemical and 
thermal equilibria with a baryon chemical potential /is = 0 (which is a good approximation 
for LHC energies) until the temperature when such an approach breaks down. Then the 
system has lost the properties of local equilibrium, thermal and chemical as well, and another 
approximation should be used. 


D. Particlization stage 


As discussed in the Introduction, the basic HKM describes particlization as a continuous 
process (for review see [l^): particles gradually “escape” from the expanding fluid, forming 


the non-equilibrium Wigner function, whic 
as the input for UrQMD hadronic cascade 
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len can be used at some space-like hypersurface 
13| . As it was already discussed, at smooth initial 
conditions the results are very similar to a sudden particlization scenario and we will use 
the latter in this paper. We assume that the chemically and thermally locally equilibrated 
evolution takes place until temperature T = 165 MeV (corresponding to an energy density 
e = 0.5 GeV/fm^ for the Laine-Schroeder EOS) is reached, and switch to particle cascade 
at the hypersurface defined by this criterion. Such a switching surface is built during the 










hydrodynamic evolution with the help of the Cornelius 1^ routine. We apply the Cooper- 


Frye formula to convert the fluid to the cascade of particles: 



( 12 ) 


In order to take into account the viscous corrections to the distribution function, Grad’s 
14-moment ansatz is used. We assume that the corrections are the same for all hadron 
species. Then formula fll2p transforms to 



This distribution function is used to create an ensemble of particles on the hypersurface. At 
hrst, the average number of hadrons of every sort is calculated: 


(14) 


ANi = Aa^u^Ui^th = 


The average total number of particles then equals (Wot) = SjW- The exact total number of 
the particles to be created is sampled according to Poisson distribution with a mean value 
< Ntot >• The type of each generated particle is chosen randomly based on probabilities 
Ni/Ntot- Then the momentum is assigned to each particle in the local rest frame of the 
fluid. The direction of momentum is chosen randomly in the 47r solid angle and its modulus 
is generated according to the isotropic part of Eq. flTdl) . After that, the corrections for 
Wresiduai ' WAisc are applied via the rejection sampling. The particle position is set to be 
equal to the position of centroid of the surface element, and the spacetime rapidity is sampled 
randomly within the longitudinal size of the volume element. Finally, the particle momentum 
is Lorentz boosted to the center-of-mass frame of the hreball. 

E. Hadronic cascade 

The generated hadrons are then fed into the UrQMD cascade. Since the cascade ac¬ 
cepts only a list of particles at an equal Cartesian time as an input, the created particles 
are propagated backward in time to when the hrst particle was created. The particles are 
not allowed to interact in the cascade until their trajectories cross the particlization hy¬ 
persurface. The Laine-Schroeder EoS, that is applied in our analysis corresponds to an 
equilibrium hadron-resonance gas consisting of about 360 hadrons in the low temperature 
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limit. Many of those heavy hadrons are not included in the UrQMD hadron list. To prevent 
violation of the energy-momentum conservation even at space-like parts of the isotherm 
T = 165 MeV, we decay the heavy resonances, which are not in the UrQMD list, just 
at the switching hypersurface. The particle propagation is stopped at Cartesian time 400 
fm/c, where their coordinates and momenta are recorded. We generate 50000 GLISSANDO 
events for each centrality class to produce averaged initial energy density prohle for the pre- 
thermal/hydrodynamic evolution. Only the one hydrodynamic calculation is performed for 
each centrality class with a given set of parameters. For each hydrodynamic evolution 20000 
UrQMD events are generated. The generated sets of events are stored in ROOT trees and 
hies, which are further processed with scripts to plot momentum distributions and calculate 
how coefficients or correlation functions for physical analysis. 


III. RESULTS AND DISCUSSION 


To perform calculations of the spectra, anisotropic how coefficients, and interferometry 
radii we need to specify the initial states and model parameters. As discussed above, we 
generate these initial conditions with GLISSANDO 2, a Monte Garlo Glauber generator. The 
output of this generator is the relative deposited strength (RDS), which dehnes the initial 
energy density ([I]) through the number of wounded nucleons Nw and number of binary 
collisions Nun- The relative weight of these quantities a is crucial for the description of all 
charged-hadron multiplicities at diherent centralities. It is found that the value a = 0.24 
at diherent reasonable values of other parameters gives the best ht for the multiplicities at 
LEG energy ^/s = 2.76 TeV. Figure [T] demonstrates a typical ht at this value of a. 

To analyze the ehects of momentum anisotropy of the initial state, similarly to Ref. Q 
we choose the momentum dependence of the distribution function fo(p) in Eq. ([3]) in the 


form 


29|: 


fo(p) = gexp 



(p-u)^-(p-vp , (p-vy 

^ o I 


A 


(15) 


where = (cosh? 7 , 0, 0, sinhr;), = (sinhr^, 0, 0, coshr;). Note, that in the rest frame of 
the huid element, rj = 0, {p ■ UY — {p ■ VY = p]_ and {p ■ VY = Then one can see 
that the parameters Ay and can be associated with the diherent temperatures along the 
beam axis and orthogonal to it correspondingly. The parameter A = A_l/A|| is the main 


10 







parameter, that defines momentum anisotropy of the initial state. We use two values for 
it: A = 1 (momentum isotropic case) and A = 100 (anisotropic case, almost no pressure in 
the longitudinal direction). We put = 1.4 GeV H- At rr^/\\ « 1 the last value is 
just absorbed into a factor g = 5 '(Ax,A||) that is the normalization constant in Eq. flTSD . 
This constant is defined from the requirement that the initial energy-momentum tensor zero 
component, formally associated with the dimensional function /o(p) only and so calculated 
according to Eq. (jl]) with substitution f{x,p) —)■ /o(p), is unity: T°°[/o(p)] = 1 [see Eqs. 
([I]) and ([3])]. As for the anisotropy factor A, it also does not affect the initial energy density 
distribution ([T]), but reveals itself in the process of the matter evolution at the pre-thermal 
stage, see Sec. II B. 

The generated distribution functions are related to the two initial times Tq = 0.1 and 
tq = 0.5 fm/c for comparison. The thermalization time is chosen to be, as discussed before, 
Tth = 1 fm/c. When other parameters are fixed, we make calculations for 0-5%, 5-10%, 
10-20%, 20-30%, 30-40%, 40-50%, and 50-60% centralities with the same normalization 
coefficient eo in Eq. ([T]) and fixed relative contribution of binary collision a = 0.24. The 
latter, as we discussed above, does not in fact depend, on all the other parameters, contrary 
to normalization eo. The eo values at different scenarios, that are associated with iHKM 
parameters, are shown in Table 1. 


model 

A 

'^rel 

v/S 

To 

(^) 

Co [GeV/fm^] 

hydro 

- 

- 

0 

0.1 

5.16 

1076.5 

hydro 

- 

- 

0.08 

0.1 

6.93 

738.8 

iHKM 

1 

0.25 

0.08 

0.1 

3.35 

799.5 

iHKM 

100 

0.25 

0.08 

0.1 

3.68 

678.8 

iHKM 

100 

0.75 

0.08 

0.1 

3.52 

616.5 

iHKM 

100 

0.25 

0.2 

0.1 

6.61 

596.9 

iHKM 

100 

0.25 

0.08 

0.5 

5.36 

126.7 


Table 1. The maximal initial energy densities eo and the mean /ndf for different scenar¬ 
ios. The values To, Trei correspond to fm/c. The mean x^/ndf is calculated by an averaging 
X^/ndf values taken from model comparison with observed pion, soft pion, kaon, and an¬ 
tiproton spectra, mid-rapidity charged multiplicity dependence on centrality, V 2 transverse 
momentum dependence for charged particles, and the pion interferometry radii as functions 
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of the transverse momentum. 

For each set of parameters with corresponding cq (Table 1), the experimental multiplicity 
dependence on centrality is reproduced well, similar to Fig. [H with the same a = 0.24. 
There is a correspondence between the ratios of eo at the different parameters and inverse 
ratios of the e nerg y densities at the thermalization time Tth = 1 fm/c, obtained at the such 


parameters in 


16l | when the initial energy densities at Tq are the same. 


After specifying the initial conditions we run the relaxation model, which describes the 
pre-thermal stage with different relaxation times r^ei = 0.25 and 0.75 fm/c. The grid spac¬ 
ing and time step for the relaxation model and further pure hydrodynamic calculations are 
Sx = 6y = 0.2 fm/c and 6t = 0.05 fm/c. Here we perform 2-|-1-dimensional longitudinally 
boost-invariant calculations in iHKM, which is a good approximation for the central rapidity 
interval at LHC energies. At the thermalization (proper) time Tth = 1-0 fm/c, the evolution¬ 
ary equations at the pre-thermal stage are smoothly switched (automatically) to relativistic 
viscous hydrodynamic equations in the Israel-Stewart framework. 

We use the two values of the shear viscosity to entropy ratio in the hydrodynamic phase: 
the minimal one y/s = 0.08 ~ ^ and y/s = 0.2 for comparison. Also we compare the basic 
scenario (BS) with tq = 0.1, Trei = 0.25 fm/c, rj/s = 0.08, anisotropy parameter A = 100, 
which is selected as giving the optimal description of the experimental data (see below), 
with results of calculations with other different parameters, including viscous and ideal pure 
hydrodynamic scenarios (without the pre-thermal stage but with the subsequent hadronic 
cascade). 

In Fig. [2] we compare the results with the pion spectra at 0.1 < pr < 3 GeV/c. As one 

can see at fairly large px > 1-5 GeV /c the results for pion transverse spectra in the scenario 

with initial time tq = 0.5 fm/c describe poorly the experimental data: xTo„/ndf=14.33. 

The best plots correspond to shear viscosity value p/s = 0.2, xTQ^/ndf=2.48, in the basic 

scenario (BS) with p/s = 0.08 —?■ 0.2. The BS itself leads to satisfactory results for pion 

2 

spectra also. According to Table 1 the basic scenario leads to mean value (^)=3.68 over 

2 

all observables under consideration, while the scenario with p/s = 0.08 —)■ 0.2 has ( /^ ) = 
6.61 because it does not describe V 2 (see below Fig. ED, x2^/ndf=27.92. 

In addition to the basic scenario the relatively good description of the wide class of 
observables gives, according to Table 1, the similar scenarios with, however, larger relaxation 
time, Trei = 0.25 fm/c —)■ 0.75 fm/c, and the scenario with isotropic initial state, A = 100 —)■ 


12 



1. So, for more detailed analysis, let us consider the interesting phenomenon of soft pion 
enhancement presented in Fig. [3l As one can see, the iHKM scenarios with isotropic initial 
state, A = 1 or with larger relaxation time, r^ez = 0.75, obviously underestimate pt spectra 
in the soft region, namely, for the isotropic case and xLftpions/’^'^f=l-27 

for Trei = 0.75 fm/c; pure (viscous) hydrodynamics scenario with Tth = Tq = 0.1 fm/c is 
also less effective at soft pt for pions. At the same time, the basic scenario leads to twice 
less value xLftpions/'^'^f~0-^9. The description of this soft momentum region becomes the 
best when the initial state formation is ascribed to the later time tq = 0.5 fm/c. However, 
it excludes, as we discussed above (see also Table 1), a satisfactory spectra description at 
Pt > 1.5 GeV/c despite that there is no problem with V2 coefficients and the interferometry 
radii. This could mean the necessity of including mini-jets in the hydrokinetic picture in 
order to describe the large px spectra simultaneously with very small transverse momenta. 
One of the other proposals considers the specific mechanism of soft pion radiation, namely, 
the Bose-Einstein condensation in the thermal model 30|. However, since the pion spectra in 
the basic scenario are within the experimental errors with small xLftpions/i^df , we consider 
now the BS as the most realistic, allowing to describe well the bulk observables in wide px 
region. 

As one can see from Fig. 01 iHKM also describes well the kaon spectra in the basic 

scenario, x^^^„/ndf=1.9. As for the antiproton spectra (Fig. [5]), one can see that the basic 

2 

scenario is not the best one (Fig. [5]). That is why, if one compares the mean values of (^) 
over all observables in Table 1, this scenario brings slightly worse results than scenarios with 
an isotropic initial state and anisotropic one with larger relaxation time. At the same time 
the last two scenarios do not describe the enhancement of soft pions as we already discussed. 

Let us apply iHKM to a description of elliptic flows. The elliptic flows in iHKM are 
associated with the reaction plane that is well-defined in GLISSANDO 2. We compare the 
results with the experimental n 2 -coefficients estimated with four-particle cumulants mea¬ 
sured for unidentified charged particles as a function of transverse momentum for various 
centrality classes, ^2(4} (311 • Figure 0] demonstrates all charged-particle n 2 -coefficients for 
the basic scenario in comparison with iHKM results with p/s = 0.2 instead oi p/s = 0.08 
and also with ideal hydro starting at Tth = tq = 0.1 fm/c. One can see that the two last 
scenarios lead to poorly described n 2 -coefficients. At the same time the basic scenario, as 
well as the iHKM one with tq = 0.5 instead of Tq = 0.1 fm/c and also the viscous pure 
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hydrodynamics scenario with Tth = tq = 0.1 fm/c describe the data quite satisfactorily, see 
Fig. [3 However, the two last scenarios describe the spectra worse: the one with tq = 0.5, 
at fairly large pt] the other, the pure viscous hydro, at both large and small momenta. 

An important observation for evolutionary models is the interferometry radii that reflect 
the space-time structure of the particle emission from an expanding hreball. The detailed 
study within iHKM shows that if a selection of some set of parameter values, such as Tq, 
Trei, A, rj/s, is accompauied by the renormalization of the maximal initial energy density cq, 
similar to that in Table 1, to keep the iHKM charged-particle multiplicity to be equal to 
the experimental one, the interferometry radii change within a few percent only, which is a 
much smaller change than the one caused by a change in the centrality class of the collision. 
Published results (see, e.g. review [l^) have shown that the interferometry radii depend on 
particle species, multiplicity, and initial system size. The other details are not so important. 

In Figs. IHHTOlwe present the pion interferometry radii in the iHKM basic scenario for three 
different centralities. In general the radii are described well except some small deviation for 
Rside and Rout from the central experimental points; these deviations are oppositely directed 
for these transverse radii, which results in slight overestimation of the Rout to Rside ratio, as 
demonstrated in Fig. [HI 


IV. SUMMARY 

The integrated hydrokinetic model (iHKM) of A+A collisions is developed. It includes 
the generation of the initial (generally momentum anisotropic) state, matter evolution at the 
pre-thermal stage leading to thermalization, subsequent viscous hydrodynamic expansion, 
particlization, and hadronic cascade UrQMD. This model is applied to describe at different 
centralities the mid-rapidity charged-particle multiplicities, the pion, kaon, and antiproton 
spectra, the V 2 coefficients for all charged hadrons, and the pion interferometry radii in Pb- 
Pb collisons at LHC with ^/snn = 2.76 TeV. It is found that a quite satisfactory description 
of the bulk observables can be reached when the initial states are attributed to the very 
small initial time 0.1 fm/c, the pre-thermal stage (thermalization process) lasts at least 
until 1 fm/c, and the shear viscosity at the hydrodynamic stage of the matter evolution 
has its minimal value, p/s = The requirement of satisfactory description of the soft 
pion enhancement phenomenon also leads to the further discrimination between the iHKM 
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scenarios, namely, the initial state should be maximally momentum anisotropic and have 
a small mean relaxation time, Trei Tth- In this basic scenario the bulk observables are 
described in a wide pt region, including the pion and antiproton spectra, as well as their 

Q 'ds when annihilation processes at the afterburner (UrQMD) stage are taken into account 

It is worthy noting that while the found tendencies, which direction the changes of spec¬ 
tra and V 2 values go when changing the basic model parameters, are probably stable, an 
accounting of event-by-event fluctuations, as well as modihcation of such important factors 
as EoS and particlization temperature, can change the numerical values of the model pa¬ 
rameters, ensuring an optimal description of the observables. Such an investigation is the 
subject of subsequent work. 

Another important point which we would like to stress is the approximate similarity 
of the results at different reasonable values of the main parameters of the model, even 
including the pure hydrodynamic scenario starting from Tth —)■ Tq = 0.1 fm/c. The reason 
for such similarity is probably that each change of basic parameters is accompanied by a 
re-normalization of the maximal energy density of the initial state, to keep the agreement 
with the experimental multiplicity of all charged particles at midrapidity. Then, at the 
same particlization temperature for all the scenarios, the main characteristics of the bulk 
observables are approximately preserved. Of course, some details are different, as iHKM 
demonstrates, but the variations are not so dramatic. This observation explains satisfactory 
agreement of different variants of hydrodynamic and hybrid models with the experimental 
data for A-l-A collisions, especially when not all bulk observables are considered. 

In this paper the Monte Carlo Glauber model is used for the generation of the initial state. 
The prohle from GLISSANDO is associated here with the one for energy density of the non- 
equilibrium and non-thermal initial state. Note that the entropy density cannot be used to 
characterize such a state at Tq = 0.1 fm/c. The interesting point is that with such treatment 
of the Monte Carlo Glauber model, the share of the binary collisions, contributed to the total 
energy density, is constant at about 1/4 for LHC energies at different iHKM parameters. We 
are planning to utilize other initial state models, such as EPOS and IP-Glasma, which have 
other structures of the initial state and already fixed values of the maximal initial energy 
density, as the generators of the initial states in iHKM. The comparison of the models have 
to be done on an event-by-event basis and include besides spectra, their anisotropy and 
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femtoscopic scales for different hadron pairs, as well as fluctuations of the observables. 


ACKNOWLEDGMENTS 

The research was carried out within the scope of the EUREA: European Ultra Relativistic 
Energies Agreement (European Research Group: “Heavy ions at ultrarelativistic energies”), 
and is supported by the National Academy of Sciences of Ukraine (Agreements F7-2016). 
I.K. acknowledges hnancial support by the Helmholtz International Center for FAIR and 
the Hessian LOEWE initiative. 


[1] R. Derradi de Souza, T. Koide, T. Kodama, Prog. Part. Nucl. Phys. 86, 35 (2016); 

[2] U. Heinz, J. Phys. Conf. Ser. 455, 012044 (2013); U. Heinz, and R. Snellings, Annn. Rev. Nncl. 
Part. Sci. 63, 123 (2013) [arXiv:1301.2826j; C. Gale, S. Jeon, B. Schenke, Int. J. Mod. Phys. A 
28, 1340011 (2013); P. Hnovinen, Int. J. Mod. Phys. E 22, 1330029 (2013) [arXiv:1311.1849j; 
H. Song, arXiv:1401.0079 [nncl-th] (2014) 

[3] M. Alvioli, H.-J. Drescher, M. Strikman, Phys. Lett. B 680, 225 (2009); H. Holopainen, H. 
Niemi, K. J. Eskola, Phys. Rev. C 83, 034901 (2011); M. Alvioli, H. Holopainen, K. J. Eskola, 
M. Strikman, Phys. Rev. C 85, 034902 (2012). 

[4] D. Kharzeev, C. Lourenco, M. Nardi, H. Satz, Z. Phys. C 74, 307 (1997); D. Kharzeev, M. 
Nardi, Phys. Lett. B 507, 121 (2001); D. Kharzeev, E. Levin, Phys. Lett. B 523, 79 (2001); 
D. Kharzeev, E. Levin, M. Nardi, Phys. Rev. C 71, 054903 (2005). 

[5] K. Werner, Phys. Rep. 232 87 (1993). 

[6] R. Paatelainen, K.J. Eskola, H. Niemi, K. Tuominen, Phys. Lett. B 74, 126 (2014). 

[7] H. Niemi, K. J. Eskola, R. Paatelainen, arXiv:1505.02677 [hep-ph] (2015). 

[8] B. Schenke, P. Tribedy, R. Venugopalan, Phys. Rev. Lett. 108, 252301 (2012); Phys. Rev. C 
86, 034908 (2012). 

[9] B. Schenke, S. Jeon, and C. Gale, Phys. Rev. Lett. 106, 042301 (2011). 

[10] S.V. Akkelin, Yu.M. Sinyukov, Phys. Rev. G 81, 064901 (2010). 

[11] Yu.M. Sinyukov, S.V. Akkelin, Y. Hama, Phys. Rev. Lett. 89, 052301 (2002). 

[12] S.V. Akkelin, Y. Hama, In.A. Karpenko, Yu.M. Sinyukov, Phys. Rev. G 78, 034906 (2008); 


16 



Yu.M. Sinyukov, S.V. Akkelin, lu.A. Karpenko, Y. Hama, Acta Phys. Polon. B 40, 1025 
(2009) [arXiv:0901.1576j; Yu.A. Karpenko, Yu.M. Sinyukov, J. Phys. G 38, 124059 (2011) 
[arXiv:1107.3745j; 

[13] lu.A. Karpenko, Yu.M. Sinyukov, K. Werner, Phys. Rev. C 87, 024914 (2013). 

[14] Yu.M. Sinyukov, S.V. Akkelin, lu.A. Karpenko, V.M. Shapoval, Advances in High Energy 
Physics. Volume 2013 (2013) [http://dx.doi.org/10.1155/2013/198928l. 

[15] P. Huovinen, H. Petersen, Eur. Phys. J. A 48, 171 (2012) [arXiv:1206.3371J; S. Pratt, Phys. 
Rev. C 89, 024910 (2014); D. Molnar, Z. Wolff, arXiv: 1404.7850. 

[16] V.Yu. Naboka, S.V. Akkelin, lu.A. Karpenko, and Yu.M. Sinyukov, Phys. Rev. C 91, 014906 
(2015), arXiv; 1411.4490 [nucl-th] (2014). 

[17] P. Romatschke and M. Strickland, Phys. Rev. D 68, 036004 (2003). 

[18] M. Strickland, Nucl. Phys. A 926, 92 (2014); arXiv: 1410.5786. 

[19] M. Martinez and M. Strickland, Phys. Rev. C 81, 024906 (2010); Nucl. Phys. A 848, 183 

( 2010 ). 

[20] W. Elorkowski, R. Ryblewski, and M. Strickland, Nucl. Phys. A 916, 249 (2013); Phys. Rev. 
C 88, 024903 (2013). 

[21] L. Tinti, W. Elorkowski, Phys. Rev. C 89, 034907 (2014); D. Bazow, U.W. Heinz, and M. 
Strickland, Phys. Rev. C 90, 054910 (2014); L. Tinti, arXiv: 1411.7268. 

[22] S.A. Bass et al, Prog. Part. Nucl. Phys. 41, 255 (1998); M. Bleicher et al, J. Phys. G 25, 
1859 (1999). 

[23] A. Krasnitz, Y. Nara and R. Venugopalan, Nucl. Phys. A727, 427 (2003). 

[24] T. Lappi, Phys.Lett. B643,ll (2006). 

[25] W. Broniowski, M. Rybczyiiski, P. Bozek, Gomput. Phys. Commun. 180, 69 (2009) 
[arXiv:0710.5731J. 

[26] M. Rybczynski, G. Stefanek, W. Broniowski, P. Bozek Gomput.Phys.Gommun. 185 (2014), 
pp. 1759-1772; arXiv: 1310.5475. 

[27] M. Laine and Y. Schroder, Phys. Rev. D 73, 085009 (2006). 

[28] lu. Karpenko, P. Huovinen, M. Bleicher, Gomput. Phys. Gommun. 185, 3016 (2014) 
[arXiv: 1312.4160]. 

[29] R. Ryblewski, W. Elorkowski, Acta Phys. Polon. B 42, 115 (2011). 

[30] V. Begun, W. Elorkowski Phys. Rev. C 91, 054909 (2015). 


17 


[31] B. Abelev et. al. (The ALICE Collaboration), Phys. Lett. B 719, 18 (2013). 

[32] K. Aamodt et al. (ALICE Collaboration), Phys. Rev. Lett. 106, 032301 (2011). 

[33] B. Abelev et al. (The ALICE Collaboration), Phys. Rev. C 88 044910 (2013). 

[34] J. Adam et al. (The Alice Collaboration), arXiv:1507.06842 [nucl-th]. 


18 



FIG. 1. Multiplicity centrality dependence for the iHKM basic scenario with the parameter values tq = 0.1 fm/c, = 0.25 
fm/c, r{/s = 0.08, and A = 100 for the LHC energy y/SNjsf = 2.76 TeV. Experimental data are from Q. 



FIG. 2. Resulting pion spectra in the 0.1 < pt < 3 GeV/c region for various centrality classes obtained in the iHKM basic 
scenario (as in Fig. The results are compared with those in iHKM at the other parameter ro = 0.5 fm/c and with pure 
viscous hydro at the starting time tq = 0.1 fm/c for centrality classes 0-5%, 10-20%, and 30—40%. The experimental 

data are from 0- The spectra for different centralities are multiplied by a factor of 2 (2® = 64 for 0—5% centrality). The 
lines correspond to different centrality classes as marked by letters: a, 0-5%; b, 5-10%; c, 10—20%; d, 20—30%, e, 30-40%; f, 
40-50%; g, 50-60% 
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FIG. 3. Detailed picture of pion spectra in soft pr region for 0-5% centrality in the iHKM basic scenario 
(as in Fig. [1} in comparison with the results obtained with (1) the other relaxation time Trei = 0.75 fm/c 
instead of 0.25 fm/c, or with (2) the isotropic parameter A = 1 instead of the anisotropy one A = 100, or 
with (3) the other initial time tq = 0.5 fm/c instead of 0.1 fm/c. Also the results for pure viscous hydro, 
starting at Tth —>■ tq = 0.1 fm/c are presented. The experimental data are from 
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[g-AQOKLtP dp diisVNjP [g-ASOKLiP dp diisVN^P 



FIG. 4. The resulting kaon spectra under the same conditions as in Fig. [2j 



FIG. 5. The resulting antiproton spectra under the same conditions as in Fig. m 
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iHKM, ys^o.os, 0-5% 
iHKM, yS^O.OS, 5-10% 
iHKM, yS^O.OS, 10-20% 
iHKM, yS^O.OS, 20-30% 
iHKM, yS^O.OS, 30-40% 
iHKM, yS^O.OS, 40-50% 
iHKM, ti/S=0.2, 5-10% 
iHKM, ti/S=0.2, 10-20% 
iHKM, ti/S= 0.2, 20-30% 
hydro, t|/S= 0, 5-10% 
hydro, t|/S= 0, 10-20% 
hydro, t|/S=0, 20-30% 


FIG. 6. All charged-particle V 2 coefBcients for centrality classes 0-5%, 5-10%, 10-20%, 20-30%, 30—40%, 40-50% and 50- 
60%, obtained in the iHKM basic scenario (as in Fig. m- The results are compared with those in iHKM at the other parameter 
dissipation condition, ry/s = 0.2 instead of 0.08 and with ideal hydro with the starting time = tq = 0.1 fm/c for centrality 
classes 5-10%, 10-20%, and 20-30%. The experimental data are from [^. 



iHKM, i:o=0.1,0-5% 
iHKM, Xo=0.1,5-10% 
iHKM, Xo=0.1, 10-20% 
iHKM, tg^O.I, 20-30% 
iHKM, Xo=0.1,30-40% 
iHKM, Xo=0.1,40-50% 
hydro, To= 0.1,5-10% 
hydro, To= 0.1, 10-20% 
hydro, To=0.1,20-30% 
iHKM, i:o=0.5, 5-10% 
iHKM, Xo=0.5, 10-20% 
iHKM, i:o=0.5, 20-30% 


FIG. 7. Same as Fig. m but the results are compared with those in iHKM at the other initial time, tq = 0.5 fm/c instead of 
0.1 fm/c and with viscous hydro at the starting time = tq = 0.1 fm/c. 
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FIG. 8. (Color online) The Rout dependence on transverse momentum for different centralities in 
the iHKM basic scenario under the same conditions as in Fig. [TJ 



FIG. 9. (Color online) The Rgide dependence on transverse momentum for different centralities in 
the iHKM basic scenario under the same conditions as in Fig. [TJ The experimental data are from 
34|. 
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FIG. 10. (Color online) The Riong dependence on transverse momentum for different centralities 
in the iHKM basic scenario under the same conditions as in Fig. [TJ The experimental data are 
from 34( 1 . 



FIG. 11. (Color online) The Rout/Rside ratio for different centralities in the iHKM basic scenario 


under the same conditions as in Fig. [TJ The experimental data are from 


m- 
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